Phonon surface mapping of graphite: disentangling quasi— degenerate phonon 

dispersions 
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The two-dimensional mapping of the phonon dispersions around the K point of graphite by 
inelastic x-ray scattering is provided. The present work resolves the longstanding issue related 
to the correct assignment of transverse and longitudinal phonon branches at K. We observe an 
almost degeneracy of the three TO, LA and LO derived phonon branches and a strong phonon 
trigonal warping. Correlation effects renormalize the Kohn anomaly of the TO mode, which exhibits 
a trigonal warping effect opposite to that of the electronic band structure. We determined the 
electron-phonon coupling constant to be 166(eV/A)^ in excellent agreement to GW calculations. 
These results are fundamental for understanding angle-resolved photoemission, double-resonance 
Raman and transport measurements of graphene based systems. 



The lattice dynamics of most "standard" materials are 
well known nowadays and can be routinely described 
by ab-initio methods based on density-functional per- 
turbation theory [1 . The phonon dispersion relations 
of graphite and graphene have already been determined 
by inelastic neutron scattering [2], inelastic x-ray scat- 
tering (IXS) [3 13], double-resonance Raman scatter- 
ing [5l E] and electron energy loss spectroscopy [7] . After 
the first ab-initio calculations [HI [3] had been corrobo- 
rated by IXS [21 SI HO], the story for graphite seemed 
closed [TT]. However, the phonon dispersions of graphite 
and graphene - a single layer of graphite - have contin- 
ued to present surprises in recent years: (i) two Kohn 
anomalies [T^] at the Brillouin zone (BZ) center (F) 
and corner (K) were predicted theoretically [T^l [S] and 
(ii) for charged graphene/graphite, the (electronic) non- 
adiabatic effects at F were predicted theoretically [TS] 
and confirmed experimentally by Raman scattering |16l 
[17] . Very similar observations have been made for car- 
bon nanotubes [HI [19j [2Ql EU [22|. The precise under- 
standing of the Kohn anomaly is highly significant for 
electron-phonon coupling (EPC) and thus of prime im- 
portance for the quantitative description of superconduc- 
tivity in graphite intercalation compounds |23l 124] , elec- 
tronic transport in the high bias regime, and for the (dou- 
ble) resonant Raman scattering [3 US] . The latter is now 
commonly used to distinguish single-layer graphene from 
double- and multi-layer graphene [^ \Tl\ . 

Qualitatively, the Kohn anomaly of the transversal 
optical (TO) branch at the K point is well described 
by density functional theory (DFT) [14|. However, its 
magnitude is severely underestimated by DFT [55] US] ■ 



Electron-electron correlation leads to an enhancement of 
the EPC at K. This can be quantitatively calculated [59] 
with the GVF-approximation which also gives very good 
results for the electronic band-structure of graphite (see, 
e.g., Ref. [3D|). Although the EPC at F can be derived 
from the IXS and Raman linewidth measurements |31j 
the experimental determination of the EPC close to the 
K point [29] is still missing. The experiments needed 
to unravel the details of the phonon dispersions are ex- 
tremely challenging because three phonon branches are 
overlapping in a small energy window and the phonon 
energies have a the strong dependence on doping. Thus, 
in order to overcome this problem there is an urgent need 
to map the full two-dimensional phonon dispersion rela- 
tions around K point in a similar manner as done rou- 
tinely for electron energy dispersions by angle-resolved 
photoemission spectroscopy (ARPES). 

In this Letter we present a mapping of the detailed 
two-dimensional phonon dispersion relation around the 
K point of graphite single crystals measured by IXS. The 
phonon dispersion perpendicular to the graphene lay- 
ers in graphite is practically zero for the phonon modes 
around K [Tl] [32] and thus we can use graphite to test 
the phonons of graphene. We prove that at K, the two 
branches with LA and LO character and the branch with 
TO character are almost triply degenerate due to the 
Kohn anomaly which brings down the TO branch in en- 
ergy. The phonon mode assignment is carried out by 
explicit comparison of the measured and calculated IXS 
intensities in the two-dimensional map and used to dis- 
entangle overlapping phonon modes. Using the present 
unambigous assignment of the phonon branches and pre- 
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FIG. 1: (a) Scattering geometry for the IXS experiments; 
spectra are collected along a line starting from the (010) 
BZ along the (210) [i.e. FKM] direction, (b) IXS spectra 
of graphite single crystals along the TKM line. The three 
phonon modes are fit by Lorentzians. The magenta and green 
Lorentzians denote the LA and LO modes. The TO mode is 
denoted by a red Lorentzian. (c) The calculated (left panel) 
and measured (right panel) IXS intensites for LO (green) and 
TO (red) modes. 



viously measured electron dispersions around K we de- 
termine the value of the EPC entirely from experiments 
and find quantitative agreement to GW calculations. 

The experiments were performed on beamline ID28 at 
the ESRF, utilizing the silicon (999) setup which yields 
monochromatized synchrotron light of 17794 eV, provid- 
ing a total energy resolution of ^--^3.3 meV [33]. The beam 
size at the sample position was 60 x 30 ^m^. As the 
spectrometer is equipped with nine analyzers, nine IXS 
spectra were recorded simultaneously, thus allowing an 
efficient coverage of the relevant portion of the phonon 
dispersion branches. Single crystalline graphite samples 
with about 1 cm diameter and 100 /im thickness were 
mounted on a goniometer. In Fig. [ija) the scattering ge- 
ometry and the coordinate system used throughout this 
work is shown. We collected IXS spectra along the line 
starting from the (010) reciprocal lattice point in the 
(210) direction. We have purposely chosen this line since 
calculations of the IXS intensities reveal the highest in- 
tensity for the TO mode in the (010) BZ. 

The phonon frequencies are the eigenvalues of 
y/Vq/M, where 2?q is the dynamical matrix at the 
wavevector q and M is the carbon mass. The DFT 
dynamical matrix D^^^ is computed using linear re- 
sponse [1 , with the exchange-correlation of Ref. f34] 
and other computational details as in Ref. 14J. The 
Kohn anomaly of the TO branch near K is entirely de- 
termined [m [29] by the contribution (to I?q) of the 
phonon self-energy between tt bands Pq (defined in Eq.l 
of Ref. [H]). In particular, the TO branch corresponds 
to the eigenvector |pq) of Pq with the eigenvalue Pq that 



is largest in modulo, pq is proportional to the square 
of the EPC between the tt bands near the Fermi level. 
In Ref. [33] it was shown that DFT underestimates \pq\ 
and that \pq\ from the more accurate GW calculations is 



larger by a factor 



1.61. Here, following Ref. [29], 



we compute the GW dynamical matrices as 



GW 



V. 



DFT 



where A is a constant independent of q. deter- 
mines the slope of the TO branch near K. A shifts 
rigidly the TO branch with respect to the other branches. 
A cannot be computed easily ab — initio (see footnote 
24 of Ref. [in]) and is fitted to reproduce the present 
measurements. Furthermore DFT slightly overestimates 
(in a systematic manner) the phonon frequencies of the 
branches not affected by the Kohn anomaly, thus all the 
frequencies are scaled by -1.3 % to simplify the compar- 
ison with measurements. 

Selected raw IXS spectra along the TKM direction 
are shown in Fig. [l|b). These spectra contain the con- 
tributions of three phonon modes. A lineshape analy- 
sis using Lorentzians yields the position and intensity 
of each mode which are assigned to the LA, LO and 
TO branches. Note that this characterization of phonon 
eigenvectors is only strictly valid for very small k away 
from r and we identify the phonon branches in the whole 
BZ according to their characterization close to F [35] . 
From symmetry considerations LO and LA modes are 
degenerate at the K point. Our results clearly show that 
LO /LA and TO are almost degenerate at and close to the 
K point. This is in striking contrast to previously deter- 
mined phonon dispersions [31 [1] and can not be explained 
by the standard DFT phonon calculations [^ . Our mode 
assignment is based on a comparison of the measured and 
calculated IXS intensities as shown in Fig. ijc) and fur- 
ther supported by a comparison to the GW calculated 
phonon frequencies as shown later. It is clear that the 
calculated and experimental IXS intensities of LO and 
TO show the same pattern and thus we can safely assign 
the TO mode as the one with higher intensity along the 
considered path in the (010) BZ. 

From the maxima of the Lorentzians a set of ~100 
data points is obtained in the QxQy plane in a region 
of ±0.6 around the K point, which is further ex- 
tended by the application of symmetry operations. In 
Fig. [2]ja,b) the raw experimental data points around the 
K point are shown. From these data points the three- 
dimensional surfaces of the phonon dispersion relations 
are interpolated using spline functions. The interpolated 
phonon maps are used for further data analysis. The 
phonon trigonal warping effect for the LA,LO and TO 
modes around the K point is illustrated from the phonon 
equi-energy contours shown in Fig. ^c) along with the 
GW calculations. Notably, the experiments and calcu- 
lations are in excellent agreement to each other which 
further supports the mode assignment performed above. 
Interestingly, the phonon trigonal warping effect of the 
TO mode is opposite to the electronic trigonal warping 
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FIG. 2: The interpolated experimental two-dimensional phonon dispersion relation for (a) the longitudinal optic and accoustic 
phonon branches and (b) the transverse optical phonon branch that has a Kohn anomaly. Dots denote the experimental data 
taken from the maxima of Lorentzian fits [see Fig. 1(b)]. (c) Equi-energy contours from the interpolated phonon surfaces for 
the three highest phonon branches around K for energies 155 meV-170 meV. o and A denote an interpolation to the raw 
experimental data in (a) and (b). Lines are GW calculations. The present color code follows the diabatic notation: branches 
with the same color have similar phonon patterns. 



effect, i.e. the slope of the TO phonon branch is higher 
in KM direction, whereas the electronic bands have a 
higher slope in KV direction [30] • For the LA and LO 
modes, the phonon trigonal warping has the same angu- 
lar dependence as for the electronic bands. The fact that 
the equi-energy contour of the TO mode is almost cir- 
cular, is also in agreement with the angular dependence 
of the kink due to EPC in the quasiparticle dispersion 
measured by ARPES, where no change of the kink po- 
sition, i.e. the energy of the coupling phonon [36^ was 
observed. Furthermore, the observed contour for the TO 
mode can explain why the Raman linewidth of the G' 
peak depends weakly on the exciting laser energy. The 
trigonal warping of electrons and the TO phonons is op- 
posite, canceling the effects of enlarging/shrinking of the 
Raman linewidths. Furthermore it is clear that the LO 
and TO branches cross along the high-symmetry KT di- 
rection but undergo anti-crossing along all the other di- 
rections. In a final Fig.[3]we show the phonon dispersions 
along the high-symmetry lines and compare it to previ- 
ous DFT calculations and our GW results. It can be 
observed that (1) the three branches are almost touch- 
ing each other at K (2) the LO and TO branches cross 
each other in KT direction and (3) the proper inclusion 
of the electron-electron correlation (through the GW ap- 




FIG. 3: Phonon dispersion along the TKM direction. 
Points are measurements extracted from the maxima of the 
Lorentzians from Fig. 1(b). Full lines are GW calculations. 
The dashed line is the TO branch from DFT calculations. 
The color code is the same as in Fig. |2] 



proach) in the calculation of the EPC is crucial to obtain 
a good agreement with the experiment. Indeed GW cal- 
culations are in almost perfect agreement with measure- 
ments while DFT calculations severely underestimate the 
phonon slope of the TO branch and thus the EPC at K 
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and overestimate the phonon energy. For precise descrip- 
tion of the experimental phonon dispersion relation, we 
provide simple fit formulas for phonon dispersions of the 
TO and LO branches. The TO phonon dispersion is given 

by 

ETo{q,0) = 2.56cos(36l)g2_62.75(72+73.07(7+149.8 (1) 

and the LO phonon dispersion is best fit by 

ELo{q,0) = -10.25gcos4(36l) - 5.32(jcos3(36')+ 
+5.47(7 cos2 (36*) + 2&AAq^ + 37.01(?- 
-25.1V cos(36l) + 151.5. 

(2) 

Here Eto{<i) and Elo{<i) denote the phonon energies in 
meV of TO and LO branch, respectively, q is the phonon 
wavevector measured in units from K point [see Fig. 
1(a)] and the angle away from the KM direction. From 
these fits we obtain the slope of the TO branch at K, 
S!^o = 73.07meVA, fr om which we can derive an experi- 
mental value for the square of the EPC constant between 
TT bands at K, {D\)f [H]. Indeed, (see Eq. (10) of [H]), 

^""^^'^ — — ' 

where {D'j^)p is defined in [29 , and all the parameters 
on the right hand side were experimentally determined 
in this work or previously by ARPES [30j. In particular, 
hw^Q = 149.8meV is the the TO phonon energy at K, 
vp = 1.05 X lO^ms-i [201 137| the Fermi velocity, M the 
carbon atom mass and oq = 2.46A the in-plane lattice 
constant. Using Eq.(|3]), we obtain (Dj^) p = 166 (eV/A)^ 
which is in excellent agreement to the theoretical GW 
value of {D']^)p = 164(eV/A)^ [29 and is to be compared 
to much lower DFT value of {D\)p = 92(eV/A)2 gg. 



Finally, concerning the temperature dependence of the 
energy ordering and the LO-TO difference at we 
performed experiments at three different temperatures 
(300K, 150K and 15K) and find no change in the ordering 
and a slight increase in the energy difference at K with 
decreasing temperature (2.2 mcV instead of 0.7 meV). 

In summary we have disentangled overlapping phonon 
branches and mapped the 2D phonon surfaces at the 
K point of graphite using a combination of IXS and 
ah — initio calculations. Excellent agreement for the 
phonon energies was found once many-body effects at 
the GW level are considered. The three highest branches 
around K point are almost triply degenerate, in stark 
contrast to previous measurements and DFT calcula- 
tions [31 13] . Furthermore the Kohn anomaly of the TO 
phonon branch at K was directly measured by IXS and 
the EPC was determined entirely from experiments. An 
opposite trigonal warping of the TO phonon branch as 
compared to the electrons and the LO phonon branch 
was found. This is important for understanding the an- 
gular dependence of the kink in the quasiparticle dis- 
persion as measured by ARPES in graphite intercalation 
compounds and graphene and the linewidth of double- 
resonant Raman scattering. 
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